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[57] ABSTRACT 

An anemometer utilizing a radionuclide counting tech- 
nique for measuring both the velocity and the direction 
of wind. A pendulum 12 consisting of a wire 13 and a 
ball 14 with a source of radiation 15 on the lower sur- 
face of the ball is positioned by the wind. Detectors 16, 
17 and 18 are located in a plane perpendicular to pendu- 
lum 13 (no wind). The detectors are located on the 
circumference of a circle and are equidistant from each 
other as well as the undisturbed (no wind) source ball 
position. The ratio of the outputs of detectors 16 and 17, 
and the ratio of the outputs of detectors 17 and 18 are 
applied to a processor 23 to select from the data stored 
in the processor the wind velocity V and the wind 
direction 4>. 

In a second embodiment of the invention (FIG. 10) the 
pendulum consists of a wire 27, a heavy ball 24, a string 
28 and a lighter ball 25 with radiation sources on the 
lower surfaces of both balls. 

5 Claims, 11 Drawing Figures 
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FIG. 4 
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COUNTING RATE IN DETECTOR 18 (COUNTS/5 MINUTES) 


U.S. Patent May 22, 1984 


Sheet 6 of 10 


4,449,400 
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RADIONUCULIDE COUNTING TECHNIQUE FOR 
MEASURING WIND VELOCITY AND DIRECTION 

ORIGIN OF THE INVENTION 5 

The invention described herein was made by an em- 
ployee of the U.S. Government and may be manufac- 
tured and used by or for the Government for govern- 
mental purposes without the payment of any royalties 
thereon or therefor. 

BACKGROUND OF THE INVENTION 

The invention relates generally to anemometers and 
more specifically concerns an anemometer that utilizes 
a radionuclide counting technique for measuring both 
the direction and velocity of wind. 

There are several types of anemometers currently in 
use for both applied research and industrial applica- 
tions. Many types of anemometers are designed to fur- 2Q 
nish only a single component of wind velocity; cup 
anemometers are an example. However, some physical 
situations warrant knowledge of the three components 
of wind velocity. Appropriate anemometers have been 
designed to furnish this information. A discussion of the 25 
advantages and the limitations of various types of ane- 
mometers can be found in Kaimal, J.C., Sensors and 
Techniques for Direct Measurement of Turbulent 
Fluxes and Profiles in the Atmospheric Surface Layer., 
Atmos. Technol., No. 7, Fall 1975, pp. 7-14. The char- 
acteristics of some of them are briefly summarized be- 
low. 

The hot wire and hot film type anemometers, which 
depend on the cooling power of the wind, are in general 
of small size and well suited for use in a wide range of 35 
atmospheric pressures. For transient flow measure- 
ments, the hot wire anemometer has no rival. Its fre- 
quency response extends to several hundred kHz. How- 
ever, calibration is found to change with ambient tem- 
perature, with contamination of the wire, and with 40 
prolonged use. The hot wire instruments are fragile and 
not very suitable for large wind direction fluctuations. 
The hot film instruments, on the other hand, are more 
rugged and also offer the advantage of substantial sig- 
nal-to-noise ratio gain due to the high resistance of the 45 
film, though their frequency response extends only to a 
few kHz. The main use of the hot film anemometers has 
been the measurement of turbulent liquids, though they 
can be used for the measurement of turbulent fluctua- 
tions in gases as well. 50 

The three-axis sonic anemometer is based on the fact 
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Remote wind velocity measurements in the atmo- 
sphere using laser doppler methods have also been re- 
ported. These methods are nonintrusive, have advan- 
tages of remote sensing in inacessible locations, and 
provide good velocity resolution. However, they re- 
quire powerful lasers and rather delicate signal process- 
ing equipment. 

Finally, the three axis propeller anemometers and the 
propellers-on-bivane anemometers are simple in design 
and moderate in cost. However, these anemometers are 
hampered by the limitations on their response, inherent 
in their structures. 

It is the primary object of this invention to provide an 
anemometer utilizing a radionuclide counting technique 
that measures both velocity and direction of wind. 

Another object of the invention is provide an ane- 
mometer that could be portable and yet not too fragile. 

A further object of the invention is to provide a small 
economical anemometer that requires low power. 

Still another object of this invention is to provide an 
anemometer that can be used for extended periods of 
time even in remote inhospitable places which are not 
easily accessible. 

Other objects and advantages of this invention will 
become apparent hereinafter. 

SUMMARY OF THE INVENTION 

The invention consists essentially of a pendulum with 
a source of radiation on the free end of the pendulum. 
Three radiation detectors are equally spaced along the 
circumference of a circle in a plane perpendicular to the 
pendulum, when the pendulum is in its neutral position 
such that all three of the detectors are the same distance 
from the radiation source when the pendulum is in its 
neutral position. Consequently, when there is wind the 
radiation source will be moved to a position indicative 
of velocity and direction of the wind. Means are pro- 
vided for obtaining the ratio of outputs of a first and 
second of the detectors and for obtaining the ratio of the 
outputs of the second and third of the detectors. The 
two ratios are processed by a processing means to ob- 
tain the velocity and the direction of the wind. 

In a first embodiment of the invention the pendulum 
is a ball attached to the free end of the pendulum with 
a source of radiation on the ball. In a second embodi- 
ment of the invention the pendulum includes a first ball 
located between the two ends of the pendulum and a 
lighter second ball attached to the free end of the pen- 
dulum with sources of radiation on both balls. 

BRIEF DESCRIPTION OF THE DRAWINGS 


that the speed of sound increases or decreases depend- 
ing on whether the sound is traveling in the direction of 
the wind or against the wind. This anemometer simulta- 
neously measures the three components of velocity and 
seems to be sensitive only to wind components along 
the acoustic paths. Its frequency response is limited only 
by path length, and its calibration remains stable for 
extended periods. However, path averaging limits use 
of this type of anemometer to heights of 4 meters and 
above. 


FIG. 1 is a schematic diagram of a first embodiment 
of the invention; 

FIG. 2 is a schematic diagram for illustrating the 
55 effect of wind on the source counter geometry of FIG. 

l; 

FIG. 3 shows the source motion relative to the three 
counters in FIG. 1; 

FIGS. 4, 5 and 6 show plots of the counting rate of 
60 detectors 1, 2 and 3 (labeled 16, 17 and 18, respectively) 
versus source angle displacement for configuration 2 in 


In the case of anemoclinometers, which are small in 
size, the wind velocity components are determined 
from pressure differences between the various sets of 
points on a metal sphere. The frequency response is 
reasonably good, limited only by the time constant of 
pressure transducers, but the probe needs to be oriented 
into the wind. 


FIG. 3; 

FIG. 7 shows plots of the ratio of the counting rates 
of detectors 1 and 2 versus the angle between the source 
65 line and detector 2; 

FIG. 8 shows plots of the ratio of the counting rates 
of detectors 2 and 3 versus the angle between the source 
line and detector 2; 
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FIG. 9 is a block diagram of the circuitry used with 
FIG. 2; 

FIG. 10 is a schematic diagram of a second embodi- 
ment of the invention; and 

FIG. 11 is a comparison of the sensitivity of the FIG. 

1 and FIG. 10 embodiments of the invention for various 
values of wind velocity. 

DETAILED DESCRIPTION OF THE 
INVENTION 

Turning now to one embodiment of the invention 
selected for illustration in the drawings, the number 11 
of FIG. 1 designates a support for a pendulum 12. Pen- 
dulum 12 consists of a flexible wire 13 with a light, 
hollow sphere 14, such as a plastic tennis ball attached is 
to its free swinging end. Ball 14 is punctured with sym- 
metrically distributed holes or otherwise appropriately 
roughened for aerodynamic stability. A radiation 
source 15 such as for example, a one-microcurie bis- 
muth-207 electron source, is deposited on the lower 20 
surface of ball 14. A detector 16, a detector 17 and a 
detector 18 are located in a plane perpendicular to pen- 
dulum 12 when the pendulum is in its neutral position. 
Detectors 16, 17 and 18 are radiation counters. The 
neutral position of pendulum 12 is the position shown in 
FIG. 1 (the position of the pendulum with no wind). 
Detectors 16, 17 and 18 are located on a circle 19, 120° 
apart, and equidistance from the radiation source 15 
when pendulum 12 is in its neutral position. The count- 
ing rate of each of the detectors is dependent on the 
solid angle subtended by the radiation counter at the 
location of the radioactive source. 

Now suppose there is a wind in a direction parallel to 
the plane containing detectors 16, 17 and 18. Ball 14, 
and consequently radiation source 15, will be displaced 
from the neutral (no wind) position. The relative count- 
ing rates in the three detectors will depend on the new 
position of radiation source 15. As shown in FIG. 2 the 
source ball 14 (B) is displaced through an angle (0,$) to 
position B' under the influence of drag D imposed on it 
by wind velocity V. The distance d/ between the dis- 
placed position of the source ball and the three detec- 
tors are given by 

d x = BF = {[A + /(l - cos#)] 2 + P + (1) 

P sin 2 e - 2ri sin 6 cos(120 - 
d 2 = BG — {[A + - cos + P + 

P sin 2 6 — 2 rl sin 6 cos 4>}* 
d 3 = BH = {[A + /(I - cos e>] 2 + P + 

P sin 2 0 — 2 rl sm 9 cos(120 + <£)}* 

where h is the height of the undisturbed source ball 
from the counter plane; 1 is the length of the suspension 
string; r is the radius of the counter circle; and F, G, and 
H are the respective locations of detectors 16, 17 and 18. 
Thus, the counting rates C/in the three detectors will be 

C 2 =Na 2 (0,<l>)G 2 (e,<t>)€2(e,<t>) (2) 

C 3 =Na 3 (6,<l>)Gi(0,<l>)€3(O,<t>) 

where N represents the number of source radiations 
(particles) emitted per second, and a/, G/, and €/ repre- 
sent the source self-absorption correction factor, the 
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solid angle correction factor, and the radiation detec- 
tion efficiency, respectively, for the ith detector. The 
source geometrical correction factor G,{0,<j>) is related 
to di by the equation where S / 

5 represents the detector surface area normal to d/. 

By calculating the relative counting rates R g in the 
three pairs of detectors for various combinations of 0 
and <f>, it is possible to develop theoretical calibration 
curves of R#(0,<f>) versus 0 and/or <J>. Such calibration 
10 curves can then be used to determine the displaced 
source ball position under the test wind field. The direc- 
tion of the source ball displacement <f> is the same as the 
wind direction. The magnitude of the source displace- 
ment angle 0 is related to the wind speed as follows: 

tan 0— (Z>/mg) (3) 

where 

D^kiCopA)** (4) 

Substituting for D in equation (3) gives 

/ \ (5) 

where 

k = CjjpA/2mg (6) 

m = source ball mass 
Cp .=* steady sphere drag coefficient 
p — density of the atmosphere 
V = wind velocity 

The value of the constant k can be easily determined 
by using appropriate values for the various parameters 
in equation (6). Thus, for C/>=0.47 (drag coefficient for 
a smooth, unperforated, spherical source ball), 

40 p= 1.23 X 10-3 g.cm-3 (Earth), p==1.32X 10- 5 g-cm- 2 

(Mars), A— 5 cm 2 , m=10 grams, g=981 cm-sec“ 2 
(Earth) and g=375 cm-sec~ 2 (Mars), k is calculated to 
be 1.47X10- 7 cm- 2 -sec 2 for Earth and 3.85x10-9 
cm~ 2 -sec 2 for Mars. Based on a calculation of k and 
45 measurement of source ball displacement angle 0, the 
wind speed is easily calculated using equation (5). As is 
evident from equation (5), the sensitivity of the instru- 
ment would be greatly dependent on the value of the 
source displacement angle 0. The source ball mass 
50 should therefore be so adjusted as to keep 0<45\ 

It should perhaps be pointed out that a value of 
Ci)=0.47 for a spherical source ball is correct only in 
the case of an unperforated ball. If the source ball has to 
be perforated or otherwise roughened for reasons of 
55 aerodynamic stability, the value of C d could increase 
significantly, depending on the size, number, and distri- 
bution of the perforations or the surface roughness pa- 
rameter. 

Although the individual counting rates in individual 
60 counters C/and the relative counting rate in any pairs of 
counters R/ycan be easily calculated from equations (2), 
it may be more expedient to use experimentally ob- 
tained data on counting rate ratio in order to avoid 
possible complications caused by slow drifts in radiation 
65 detector efficiencies as well as changes in electronic 
system gains. Such calibration curves were developed 
using experimentally obtained counting rate data in the 
individual detectors which are disclosed in FIGS. 3-8. 
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To develop the calibration curves three Geiger- 
Miiller tubes (detectors 16, 17 and 18) were located at 
equal distances along the circumference of circle 19 
having a diameter of 0.3-m. A one-microcurie bismuth- 
207 electron source 15 was attached to the bottom of a 
0.025-m-diameter, plastic, hollow sphere ball 14. (1 
microcurie equals 3.7 XlO 4 becquerels.) The ball was 
suspended from a rod at a height of 0.15 m from the 
Geiger-Muller counter plane. The source ball was 
moved along straight lines inclined at different angles 
with respect to the three counters in order to simulate 
different wind velocities. These lines 1 through 8, called 
configurations, are shown in FIG. 3. 

Counting rates in the three detectors were measured 
for a preselected time for different source-counter con- 
figurations. FIGS. 4, 5 and 6 show individual counting 
rates in the three detectors as a function of angle of 
displacement of the source ball from the local vertical 
for configuration 2 in FIG. 3. Note that the counting 
rates for the three detectors are different at 0=0 where 
one would expect the rates to be the same. Different 
counting rates in the three detectors arise from the fact 
that they did not have identical radiation detection 
efficiencies. However, these different counting rates for 
the detectors will not affect the accuracy, since as dis- 
cussed below counting rate ratios are used. From the 
data of the type shown in these figures relative counting 
rates were determined for the detector pairs (16,17), 
(17,18), and (18,16). FIGS. 7 and 8 show relative count- 
ing rates for detector pairs (16,17) (R 12 ) and for detector 
pairs (17,18) (R 23 ) for various configurations. 

It is clear from FIGS. 7 and 8 that the source-counter 
configuration, and hence the direction of motion of the 
source ball with respect to a predetermined reference 
line, can be easily determined from the measured values 
of the relative counting rates in the three detectors. For 
example, if the ratio of the counting rates for detectors 
16 and 17 (R 12 ) is 4.18 and the corresponding value for 
detectors 17 and 18 (R 23 ) is 2.00, then the source is 
displaced through an angle 0=10° in configuration 
1(— ). Similarly, if .R. 12 — 2.00 and R 23 = 1.40, the source 
is displaced through an angle 0=2° in configuration 
2(— ). A noteworthy feature of FIGS. 7 and 8 is that 
appropriate combinations of detector channels permit a 
reasonably uniform sensitivity across the entire angular 
range <f>. 

In view of the discussion above, the outputs of detec- 
tors 16 and 17 are applied to a divider 21 to obtain the 
ratio R 12 of the output of detector 16 to the output of 
detector 17; and the outputs of detectors 17 and 18 are 
applied to a divider 22 to obtain the ratio R 23 of the 
output detector 17 to the output of detector 18. The 
ratios R 12 and R 23 from dividers 21 and 22 are applied to 
a processor 23 which stores the calibration information 
shown in FIGS. 7 and 8. The values for 0 and <f> are 
taken from FIGS. 7 and 8 for each R 12 , R 23 combination 
and the values for 0 are used to calculate corresponding 
values of V by means of equation (5). Then the values 
for V and <j> are stored in processor 23. As stated above, 
this information can also be calculated and stored in 
processor 23. In response to the measured R 12 and R 23 
values, the corresponding V and <f> values are selected 
and applied to the outputs of the processor. Processor 
23 can be a programmable read only memory (PROM) 
which is well known and commercially available. 

Alternatively, processor 23 can be a computer, which 
in response to R 12 selects the corresponding set of (0,<J>) 
values from the FIG. 7 data and in response to R 23 


selects the corresponding set of (0,<f>) values from the 
FIG. 8 data. Then the processor obtains the intersection 
of the two selected sets to provide the (0,<f>) correspond- 
ing to the R 12 and R 23 inputs. Processor 23 then calcu- 
5 lates the value of V from 0 by use of equation (5) to 
produce the values of V and <j> at its two outputs corre- 
sponding to the inputs R 12 and R 23 . A third divider 31 
can be used to provide R 31 ; the ratio of the output of 
detector 18 to the output of detector 16. R 31 may be 
10 needed to provide a uniform sensitivity across the entire 
angular range <j>. 

Certain atmospheric environments experience large 
changes in wind speed. In order to measure wind speeds 
with reasonably uniform sensitivity over these wide 
I 5 ranges with the FIG. 1 embodiment of the invention, it 
is necessary to use either a series of instruments with 
overlapping ranges or a single instrument with an ad- 
justable source ball weight and/or size. Neither of these 
two alternatives appears to be satisfactory for practical 
20 reasons. 

For use in large changes in wind speed a second 
embodiment of the invention as shown in FIG. 10 in- 
volves the use of a double pendulum incorporating a 
heavier, metallic ball 24 of the same dimensions as the 
25 lighter, plastic source ball 25. For lower velocities, the 
lighter source ball 25 will be displaced a considerable 
amount, while the heavier, metallic ball 24 will barely 
move from its normal equilibrium position. For higher 
velocities, the source ball will be displaced through 
30 almost 90% but the displacement of the heavier ball will 
be much less, thus offsetting the loss of sensitivity asso- 
ciated with larger angles of source ball deflection. By a 
proper system design, it is possible to keep the impact of 
the wake of the upper heavier ball on the lighter source 
35 ball to a minimum, thus ensuring its stability. Such a 
double pendulum system is described below. For the 
sake of simplicity, the problem of source ball displace- 
ment is discussed in one dimension only, but the treat- 
ment can be easily extended to three dimensions. 

Consider two balls of masses mi and m 2 . The first 
(heavier, metallic ball 24) is attached to a fixed support 
26 by a light inextensible string 27 of length 1|. The 
second (lighter, plastic source ball 25) is attached to the 
bottom of the first ball by a similar string 28 of length I 2 . 
45 If the two balls have the same diameter and similar 
surface roughness characteristics, the drag force on 
them will be equal, and the equilibrium conditions can 
be described by the following equations: 


40 


50 


55 


tan 


2D 


(mi + m 2 ) g 


tan 62 — 


m 2& 


where 0i is the angle of displacement of the first, 
heavier ball, and 02 is the angle of displacement of the 
second, lighter ball. 

These equations can be used to design a system that is 
60 appropriate for a wide range of wind velocities. The 
sensitivity of such a system will depend strongly on the 
rate of change of the source ball position with wind 
speed. For the one-dimensional case, one obtains 

65 (dlx/dV)—4kj\\V co£ 0i-f 2/t 2 l 2 Kcos 3 0 2 

where \ x is the displacement of the source ball from its 
normal undisturbed position and 
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1 2(m\ -f m 2 )g 


k 2 = 


2/«2g 
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In the case of a single pendulum anemometer, the corre- 
sponding equation is 

, 10 
(dlx/dV)*s2k\V cos 3 0 

FIG. 11 shows a graph of dl*/dV as a function of 
velocity for single and double pendulum arrangements 
for a specific case. It is apparent that the use of a heavy 
compensating ball considerably improved the unifor- 
mity of response of the instrument. 

It should be emphasized, however, that a part of the 
suspension string of the lower source ball may be ex- 
posed to the wake of the first ball at higher wind speeds, 
thus making the position of the source ball unstable. If 2 
proper mechanical design does not alleviate the prob- 
lem, one of the following two approaches might pro- 
vide the solution: 

(1) Use alpha counting for lower wind speeds and 
electron counting for higher wind speeds. This can be 2 
accomplished by depositing americium-241, an alpha 
source 29, on the lower (lighter) ball and bismuth-207, 
an electron source 30, on the upper (heavier) ball and 
using a retractable, thin (< 1 mg/cm 2 ) metal absorber in 
front of each detector. Both electrons and alpha parti- 30 
eles would be counted at lower speeds when the ab- 
sorber foil would be retracted, whereas only electrons 
would be counted at higher speeds when the absorber 
foil would be introduced in front of each detector. 
Thus, any instability in the position of the lower (ligh- 35 
ter), americium-241 bearing ball at higher speeds will 
not affect the counting rates, which will be affected 
only by the position of the upper (heavier) ball bearing 
bismuth-207, an electron emitter. 

(2) Use a reversible double pendulum system wherein 40 
each ball will have bismuth-207 sources of equal 
strength deposited on it. For lower speeds (5 to 20 
m/sec), the lighter ball will be suspended below the 
heavier ball, whereas these positions will be reversed at 
higher speeds (20 to 100 m/sec). By a suitable adjust- 45 
ment of the source weights and suspension string 
lengths, the lower ball can almost always be kept out of 
the wake of the upper ball. 

The invention is a technique for measuring wind 
speeds over a wide range (up to 100 m/sec). The tech- 50 
nique is based on the inverse-square-law variation of the 
counting rate as the radioactive source-to-counter dis- 
tance is changed by the movement of the source due to 
wind velocity. An anemometer based on this principle is 
usable in all noncorrosive media and in all locations, 55 
except places with no local gravitational force. It is 
suitable for measuring wind fluctuation rates of up to 1 
kHz, being limited only by the inherently slower me- 
chanical response functions of the suspension system. 

Even though the discussions presented have assumed 60 
that the radiation counters are in a plane normal to the 
local vertical and that the wind direction is parallel to 
the counter plane, cases when the counter plane is not 
normal to the local vertical and when the wind is mov- 
ing at an angle other that 90° to the local vertical can be 65 
easily treated by appropriately modifying the calibra- 
tion constant of the instrument. The presence of a verti- 
cal wind component introduces lift in addition to the 
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drag caused by the horizontal component, thus requir- 
ing a change in the calibration constant. This type of 
wind motion, however, is not usually encountered near 
terrestrial or planetary surfaces. 

Geiger-Miiller tubes and a weak electron source were 
used in the system described herein. It may, however, 
be preferable to use semiconductor counters and stron- 
ger radioactive sources for higher accuracy as well as 
higher frequency response. It may also be necessary to 
use gamma radiation sources in certain terrestrial appli- 
cations. 

What is claimed as new and desired to be secured by 
Letters Patent of the United States is: 

1. A method for measuring wind velocity and direc- 
tion comprising the steps of: 

providing a radiation source on the end of a pendu- 
lum with the pendulum being substantially perpen- 
dicular to the wind direction when in its no wind 
position and with the pendulum moving the radia- 
tion source to a position indicative of both the 
velocity and direction of the wind when there is 
wind; 

measuring the radiation from the radiation source at 
three equally spaced locations on a circle in a plane 
perpendicular to said pendulum in its stable posi- 
tion; 

obtaining two ratios of the three measurements of the 
radiation from the radiation source; and 

determining from the two ratios the velocity and 
direction of the wind. 

2. An anemometer for measuring both the velocity 
and direction of wind comprising: 

a pendulum which in its no wind position is essen- 
tially perpendicular to wind direction; 

a source of radiation on the free end of said pendu- 
lum; 

three radiation detectors equally spaced in a plane 
perpendicular to said pendulum when the pendu- 
lum is in its no wind position such that all three of 
the detectors are the same distance from the radia- 
tion source when the pendulum is in its no wind 
position for detecting the radiation from said 
source of radiation; 

means connected to the outputs of a first and second 
of said three radiation detectors for obtaining a 
ratio of outputs of said first and second detectors; 

means connected to the outputs of the second and 
third of said three radiation detectors for obtaining 
a ratio of the outputs of the second and third detec- 
tors; and 

processing means receiving said ratio of the outputs 
of the first and second detectors and said ratio of 
the outputs of the second and third detectors for 
determining the velocity and the direction of the 
wind. 

3. An anemometer according to claim 2 wherein said 
processing means includes means for storing the plot of 
ratios of the outputs of the first and second detectors 
versus wind direction and versus wind velocity and the 
plot of the ratios of the outputs of the second and third 
detectors versus wind direction and versus wind veloc- 
ity and means responsive to the ratio of the outputs of 
the first and second detectors and the ratio of the out- 
puts of the second and third detectors for selecting from 
said storing means the corresponding wind direction 
and wind velocity. 
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4 . An anemometer according to claim 2 wherein said 
pendulum includes a ball attached to its free end with 
said source of radiation on said ball. 

5. An anemometer according to claim 2 wherein said 
pendulum includes a first ball located between the two 5 


10 

ends of the pendulum and a lighter second ball attached 
to the free end of the pendulum with said source of 
radiation on said second ball. 
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